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ON THE MEAN FIELD AND CLASSICAL LIMITS 
OF QUANTUM MECHANICS 

FRANgOIS GOLSE, CLEMENT MOUHOT, AND THIERRY PAUL 


Abstract. The main result in this paper is a new inequality bearing on solu¬ 
tions of the A-body linear Schrodinger equation and of the mean field Hartree 
equation. This inequality implies that the mean field limit of the quantum 
mechanics of N identical particles is uniform in the classical limit and pro¬ 
vides a quantitative estimate of the quality of the approximation. This result 
applies to the case of C 1,1 interaction potentials. The quantity measuring the 
approximation of the A-body quantum dynamics by its mean field limit is 
analogous to the Monge-Kantorovich (or Wasserstein) distance with exponent 
2. The inequality satisfied by this quantity is reminiscent of the work of Do- 
brushin on the mean field limit in classical mechanics [Func. Anal. Appl. 13 
(1979), 115-123]. Our approach of this problem is based on a direct analysis 
of the A-particle Liouville equation, and avoids using techniques based on the 
BBGKY hierarchy or on second quantization. 


In memory of Louis Boutet de Monvel (191^1-2014) 


1. Statement of the problem 


In nonrelativistic quantum mechanics, the dynamics of N identical particles of 
mass m in R d is described by the linear Schrodinger equation 

h 2 N N 

m* = X + \ X v ( x k- x i)^ ) 

k —1 k,l =1 

where the unknown is 4' = ^(t, x \,..., xn) £ C, the A-particle wave function, 
while # 1 , 2 : 2 ,... ,Xn designate the positions of the 1st, 2nd,..., iVth particle. The 
interaction between the fcth and Itii particles is given by the potential V, a real¬ 
valued measurable function defined a.e. on R d , such that 

(1) V(z) = V(-z), for a.e. z <E R d . 

Denoting the macroscopic length scale by L > 0, we define a time scale T > 0 such 
that the total interaction energy of the typical particle with the N— 1 other particles 
is of the order of m(L/T) 2 . With the dimensionless space and time variables defined 
as 

x := x/L , t := t/T , 
the interaction potential is scaled as 


V(z) := 


NT 2 
mL 2 


V(z). 
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In terms of the dimensionless parameter 

e := HT/mL 2 

and the new unknown 

^(t,x i,... ,xn) ■= ^(t,x i,.. ,,xi v), 
the Schrodinger equation becomes 

N i N 

idi^ = + 2Ne 5Z ^(**-*0®- 

/c—1 A:,Z=1 

In the present paper, we obtain a new estimate for solutions of this Schrodinger 
equation which is of particular interest in the asymptotic regime where 

tCl (classical limit) and N 1 (mean field limit). 


Henceforth we drop all hats on rescaled quantities and consider the Cauchy problem 


( 2 ) 


N 


N 


idt^e,N = A Xk ^e,N + ^ V(x k - Xl)^ e , N , 


'I' 


«.JV|*=o 


= 'b! 


fc=i 


,N ■ 


fe,i=l 


While the discussion above applies to all types of particles, until the end of this 
section we restrict our attention to the case of bosons, i.e. to the case where the 
wave function is a symmetric function of the space variables x \,..., Xn- 
A typical example of relevant initial data for (2) is 


K n N {x 1 ,...,x N ) 


N 

( 7re ')~ dN / 4 g -(* j - 9 ) 2 / 2 « e * P -^/« ^ 

i=i 


where p,q £ R d are parameters. In other words, in this example, the initial wave 
function is the IV-fold tensor product of Gaussian wave functions with width ^/e 
and oscillations at frequency 0(l/e). 


1.1. The mean field limit. The mean field limit is the asymptotic regime where 
N —> oo, with e > 0 fixed. Set the initial data in (2) to be 


N 

*:>(*!, • • • ,x N ) := ipr(x k ) with 

fc=l 



W? (x)\ 2 dx = 1 , 


and let 'b £i jv be the solution of the Cauchy problem (2) — which exists for all times 
provided that V is such that 

N i N 

+ ^ v(x k ~ Xl ) 

fc=1 k,l =1 


has a self-adjoint extension as an unbounded operator on L 2 ((H d ) N ). Under various 
assumptions on ip m and V, it is known that, for each t £ R, 



4Uiv(f, x,z 2 ,..., 0jv)^ , e ,jv(C y,z 2 ,..., z N )dz 2 .. ■ dz N -> ip e (t, x)ip e (t, y) 
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in some appropriate sense as N —> oo, where ip e is the solution of the Hartree 
equation 

( id t %l) € = -|eA x V’e + -ipe(t,x) [ V(x - z)\ip e (t, z)\ 2 dz, 

(3) < e Jr* 

Ul,.o = <fr- 

See [29, 4, 8, 1, 9, 10, 12, 28, 26, 27, 17] for various results in this direction, obtained 
under different assumptions on the regularity of the interaction potential V. Most 
of the physically relevant particle interactions, especially the case where V is the 
Coulomb potential, are covered by these results, but not always with a quantitative 
error estimate. 


1.2. The classical limit. The classical limit is the asymptotic regime where e —> 0 
while N is kept fixed in (2) — or simply e —> 0 in (3). The formalism of the 
Wigner transform is perhaps the most convenient way to formulate this limit. Given 
<1> = <1>(X) £ C, an element of L 2 (R”), its Wigner transform at scale e is 


S) := ^ f $ (X + \eY) <l> {X - \eY)e~^ Y dY . 

J R" 

Assume that the initial data in (2) is a family 'k* rl iv such that 

We[K n N \ F$ in <S'((R d x R*)*) as e —>■ 0. 

Then, for all t £ R, the family of solutions ^ e ,N of the Cauchy problem (2) satisfies 

W e [* £)JV (i, •)] -> F N (t, •, •) in S'((R d x R d ) N ) as e -A 0 , 

where Fn = Fjv(t, Xi ,..., Xn, £i, • • •, £iv) > 0 is the solution of the following 
Cauchy problem for the IV-body Liouville equation of classical mechanics 

{ N j N 

d t F N + ' ^ x k F N --jy X! W(x k - Xi) • \7$ k F N = 0 , 

k= 1 V k,l =1 

tt * I _ cAn 

* N \t=0 “ ■ 

The classical limit of the Hartree equation (3) can be formulated similarly. Assume 
that the initial data in (3) is a family ip™ such that 

W e [ip™] -A f in in S'(R d x R d ) as e -A 0. 

Then, for all t £ R, the family of solutions ip e of the Hartree equation (3) satisfies 

W e [ip e (t, •)] -A f{t , •, •) in S'{ R d x R d ) as e -A 0 , 


where / = f{t,x,f) > 0 is the solution of the following Cauchy problem for the 
Vlasov equation of classical mechanics with interaction potential V: 


I »/ + f - V,/ - W(* - z)f(t,z)iz 'j - V s / = 0 , 

1 /l,.o = r . 

See [20, 13] for results on the classical limit of quantum mechanics involving the 
Wigner transform. 
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1.3. The mean field limit in classical mechanics. There is also a notion of 
mean field limit in classical mechanics, which can be formulated as follows. Assume 
that the initial data in (4) is 

N 

Fn(x 1, • • -,XN,€u ■ • • >&v) = f tn ( x k>€k) , 

k =1 

where f m is a probability density on R d x R d . Under various assumptions on the 
potential V, the solution Fjy of (4) satisfies 

/ F N (t,x,x 2 , ■ ■ .,x n ,£,&, ■ ■ -,£N)dx 2 d& ■ ..dx N d£ N -> f(t,x,£) 

R d xR d ) JV - 1 

in some appropriate sense as N — > oo, where f is the solution of the Cauchy problem 
for the Vlasov equation (5). See [24, 6 , 7] for the missing details. All these references 
address the problem of the mean field limit in terms of the empirical measure of the 
./V-particle system. Typically these results cover the case where V 6 C 1,1 (R), but 
the case of a Coulomb, or Newtonian interaction remains open at the time of this 
writing. For a formulation of the same results in terms of the BBGKY hierarchy, 
see [15, 22], 

Yet, there has been some recent progress on the case of singular potentials. In 
[16], the mean field limit has been established for interaction potentials with a sin¬ 
gularity at the origin that is weaker than that of the Coulomb potential. Another 
approach to the mean field limit in the case of singular interaction involves a trun¬ 
cated variant of the potential with a cutoff parameter 77 = rj(N) > 0 assumed to 
vanish as the number of particles N —> 00 : see [16, 18, 19] for the most recent 
results in that direction. This cutoff parameter can be thought of as being of the 
order of the size of the interacting particles, as explained in [18]. 

The situation described above can be summarized in the following diagram: the 
horizontal arrows correspond to the mean field limit, while the vertical arrows 
correspond to the classical limit. 


Schrodinger 


N —>-oo 


Hartree 


I 


4- 


e-S> 0 


e —>■ 0 


i 


i 


Liouville 


Vlasov 


However, these various limits are established by very different methods. The 
classical limits of either the IV-body Schrodinger equation or of the Hartree equation 
are obtained by a compactness argument and the uniqueness of the solution of 
the Cauchy problems (4) or (5). Error estimates for these limits require rather 
stringent assumptions on the regularity of the potential V and on the type of 
initial wave or distribution functions considered. The mean field limit in quantum 
mechanics (the upper horizontal arrow) comes from trace norm estimates on the 
infinite hierarchy of equations obtained from the BBGKY hierarchy in the large 
N limit. The trace norm is the quantum analogue of the total variation norm on 
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the probability measures appearing in the classical setting. In general, the total 
variation norm is not convenient in the context of the mean field limit, since it does 
not capture the distance between neighboring point particles. This suggests that 
the trace norm is not appropriate to obtain controls on the large N (mean field) 
limit which remain uniform in the vanishing e (classical) limit. 

Another notable difficulty with this problem is that the mean field limit in clas¬ 
sical mechanics is obtained by proving the weak convergence of the ./V-particle 
empirical measure in the single-particle phase space to the solution of the Vlasov 
equation. Since there does not seem to be any natural analogue of the notion 
of empirical measure for a quantum V-particle system, the analogy between the 
quantum and the classical mean field limits is far from obvious. 

Our main result, stated as Theorem 2.4 below, is a new quantitative estimate 
for the mean field limit N —> oo of quantum mechanics which is uniform in the 
classical limit e —>■ 0. 


2. Main result 

Let d be a positive integer. Henceforth we set 9) := L 2 { R d ), and 9)n '■= — 

L 2 (( R d ) N ) for each N > 1. We designate by C{$)) the algebra of bounded linear 
operators on the Hilbert space S). We denote by V(S)n) the set of operators A £ 
C(9jn) such that 

A = A* > 0 , and trace(A) = 1. 

We are concerned with the iV-body Schrodinger equation written in terms of 
density matrices, i.e. the von Neumann equation 


II 

eg 

N i AT 

2 e J2 Ak + 2 Ne ^ Vk 'AA,N 

k— 1 /c,Z=l 

l Pe,N t=0 - 

in 

Pe,N 5 


where p £> jv(f) £ V(9jn), while 

A k := ® A ® lf N ~ k) , 

and 


(7) {V jk ip)(xi,...,x N ) := V(x k - Xj)il>(xi,...,x N ) , for each ip £ Sj N . 

The notation I% designates the identity on the Hilbert space iy 

We shall everywhere restrict our attention to symmetric density matrices, cor¬ 
responding to indistinguishable particles. In other words, we assume that 

(8) T a p e , N (t)T* = Pe,N(t) 


for each t £ R and each cr £ 6^, where T a is the unitary operator defined on 
L 2 ((R, d ) N ) by the formula 


(9) T„$(X! ,...,X N ) := — 1 (1 ), . . . , 2 ^- 1 ^)) . 

One easily checks that the condition 


(10) 


T <rP<i,N T a ~ Pe,N 
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implies that (8) holds for each t G R, since 
N . N 


r ^-3 e E A * 


i 


k =1 


2 Ne 


E 


k,l=1 


= 0 for each cr € 6 


n ■ 


On the other hand, we consider the corresponding mean field equation, i.e. the 
Hartree equation written in terms of the density matrix p e (t) £ 


II 

os 

•<s> 

+ -V Pel p e 

II 

O 

p?, 


( 11 ) 

where V Pe designates both the function 

V Pe (t,x):= V(x - z)p e (t,z,z)dz 

J R d 

and the time-dependent multiplication operator defined on Sj by 


C V Pe tp){t,x ) := V Pe (t,x)ip{x). 

Next we formulate the mean field limit in terms of density operators. For each 
TV-particle density operator p n G V(S)n), we define its first n-particle marginal 
density operator, denoted by p ^ for each integer n such that 1 < n < N, by the 
following conditions: 

J G > and 

| trac e^ n {Ap'f f ) = trace i5iv ((A ® I^ N _ n )p N ) for each A G £(f) n ) ■ 

In the mean field limit, i.e. for N —> oo while e > 0 is kept fixed, one expects 
that the sequence p\ N of first marginals of the density operators p £> n solutions 
of ( 6 ) converges in some topology to the solution p £ of ( 11 ), provided that p l f n N 
approaches (p l e n )® N in some appropriate sense. The difference between p\ N and p e 
is measured in terms of a quantity analogous to the Monge-Kantorovich distance 
used in the context of optimal transport. 

First we define the notion of coupling between two density operators. 


Definition 2.1. Let d be a positive integer and let fj := L 2 (R d ). For each p, p G 
V(Sj), let <2(p,p) be the set of R G D(Sj 2 ) such that 

{ trace^ 2 (( £ 4 )R) = trac esj(Ap) 

trace^ 2 ((/f, ® A)R) = trace^(Ap) 

for each A G £(ij). 


Next we define two unbounded operators on fy 2 — T 2 (R d x R d ), as follows: 


( 12 ) 


{Qif)(xi,x 2 ) := (*1 - x 2 )if(x 1 ,x 2 ), 
{Pf)){xi,x 2 ) := -itfVu - \/ X2 )ip(xi,x 2 ), 


so that 


(13) {P*Pif)(x i - x 2 ) = —e 2 (div Xl - div X 2 )(V a;i - V X2 )ip(x 1 ,x 2 ) ■ 


The quantum analogue of the Monge-Kantorovich distance with exponent 2 is 
defined as follows. For the definition of Monge-Kantorovich distances, also called 
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Wasserstein distances, see formula (20) in section 3, or chapter 7 in [31], or chapter 
6 in [32], 

Definition 2.2. For each pdp G T>(Sj), we set 

MK 2 (p,~p) := inf tracer {{Q*Q + P*P)R) 1 ^ 2 
ReQ(p,p) 

with the following convention: 

trac e^{{Q*Q + P*P)R) := tracer (R 1/2 (Q*Q + P*P)R 1/2 ) 
if R 1 / 2 (Q*Q + P*P)R 1 / 2 is a trace-class operator, and 
trac efi 2 ((Q*Q + P*P)R) := +oo 

otherwise. 

The quantity MK 2 is not a distance on In fact, for each density operator 

p on fj, one has MK 2 (p,p) > 0 (see formula (14) below). However, MK 2 can be 
compared with the Monge-Kantorovich distance with exponent 2 (see formula (20) 
below with p = 2 for a definition of this distance, and Theorem 2.3 (2) below for 
more information on this comparison), at least for a certain class of operators and 
in the limit as e —> 0. 

Henceforth, we denote by V{R d ) the set of Borel probability measures on R d . 
For each q > 0, we define 

V q (R d ) := |/i G P(R d ) s.t. j \x\ q p(dx) < oo 

Theorem 2.3 (Properties of MK 2 ). Let d be a positive integer and let Sj := 
L 2 (R d ). For each p,p G T>(Sj) and each e > 0, one has 

(14) MK%(p,p) 2 > 2de. 

(1) Let e > 0 and let p\ and p\ be Toplitz operators at scale e on L 2 ( R d ) with 
symbols (2ne) d p\ and ( 2ire) d p2 . where pi,p 2 G V 2 (R 2d ). Then 

MK \{p\-,P 2 ) 2 < inf trace^0jj(((5*Q + P*P) OPf ((2ire) 2d p)) 

= dist M K,2(Mi,At2) 2 + 2 de. 

(2) Let p\ and p\ G 'D[Sf) > with Husimi transforms at scale e denoted respectively 
We[pf] and W t [p1f\. Then 

MK* 2 (p\,pl) 2 > dist M K, 2 {We[p\],W e [pVf) 2 - 2 de. 

Assume further that the Wigner transforms at scale e of p\ and p\, denoted re¬ 
spectively Wijpf] and Welpffl, converge in <S'(R 2d ) to Wigner measures denoted 
respectively W\ and W 2 as e —> 0. Then 

dist M K,2(wi,«;2) < lim MK 2 (p\, p 2 ) ■ 

e —^0 

The definitions and basic properties of Toplitz operators, Wigner and Husimi 
functions are recalled in Appendix B. 

Statement (2) in Theorem 2.3 implies in particular that the quantity MK 2 is not 
vanishing for all density matrices as e —> 0 + . This property is obviously essential; 
otherwise, the quantity MK 2 would not be of much practical interest for controlling 
the error in the mean field limit. Statement (1) in Theorem 2.3 will be used in 
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choosing the initial quantum states to which our error estimate for the mean field 
limit will apply. 

The main result in this paper is the following theorem. 


Theorem 2.4. Let d be a positive integer. For each e > 0 and each integer N > 1, 
let p l J) N G V(L 2 ((H d ) N )) satisfy (10) and let p™ G V(L 2 (H d )). Let t i-A p e ,N(t) be 
the solution of the quantum N-body Cauchy problem (6), and let t H > p e (t) be the 
solution of the quantum mean field Cauchy problem (11). Then, for each positive 
integer n < N and all t>0, one has 
(15) 


0 A t 


-MKftp&y 

n 


,pInV)Y <^liv^||i. 


-1 


0 A t 


A 


+ ^MKI((pTf N , 


in \2 

Pe,N> 


where 

A := 3 + 4Lip(W) 2 . 

In the particular case where p l e n N is a Toplitz operator at scale e with symbol 
(2t re) dN piy while p *" is a Toplitz operator at scale e with symbol (2 tt e) d /K™, for 
each positive integer n < N and all t > 0, one has 


(16) 


MK!)(p e {tT n ,Pl N {i)? < [2de+ A d istMK,2((Mr) 0Ar ,<iv) 2 ) e At 


8 n„ , l0 


„A t 


- 1 


A 


In particular, if p' e n G 2?(i5) is a Toplitz operator at scale e and p™ N = (p“)' s ’ iv , 
then 

(17) ^MKI(p e (tr n ,pl N (t)) 2 < (2de+ A|| W || 

Observe that one cannot deduce the mean field limit of the quantum A^-body prob¬ 
lem (6) from the bound on MK!)(p e (t ), p\ N (t)) obtained in Theorem 2.4 in the case 
where e > 0 is kept fixed, because of the term 2deexp(Af) on the right hand side 
of (15). 

On the other hand, the mean field limit alone, i.e. for e fixed, has been proved 
by other methods in this case (see [29, 4]). Moreover, quantitative estimates for 
that limit for e fixed and N —> oo have been obtained in [28, 26, 2]. Therefore, only 
the case where both N —> oo and e —> 0 remains to be treated, and the present 
work answers precisely this question. 

Indeed, the estimate in Theorem 2.4, together with the first lower bound in 
Theorem 2.3 (2), implies that the mean field limit, i.e. the convergence 

P™,N(t) Pe(t) 0n 



for each n > 1 as N —> oo is uniform as e —> 0 and over long times intervals, in the 
following sense. Let c, d satisfy 0 < c < d < 1, and set 


T(c, p, N) := — min ( In —, In N 


i(d,p,N) := 




~ 1 ,N 1 ~ C ' 


(where [x] designates the largest integer less than or equal to x ). Then the quadratic 
Monge-Kantorovich distance between the Husirni transforms at scale e of p^^(t) 
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and pf n (t) satisfies 

SUp SUp SUp dist M K,2(We[pe{t)^ n ],W e [pl N {t)]) -> 0 

1 <n<n(c',v,N) 0<t<T(c,r;,IV) 0<e<r? 

in the limit as + 77 —> 0 . 

Earlier works have discussed the mean field limit of the quantum iV-body prob¬ 
lem in the small e regime, more precisely, in the case where e = e(N) —> 0 as 
N —> 00. 

The case e(N) = V -1 / 3 is of considerable importance for the mean field limit 
of systems of N fermions, and has been investigated in [23, 30]. In the more 
recent reference [5], the iV-body problem in the fermionic case is compared to 
the Hartree-Fock equations, with error estimates in the Hilbert-Schmidt and trace 
norms — see formulas (2.19)-(2.20) in Theorem 2.1 of [5]. Of course, convergence 
in either Hilbert-Schmidt or trace norm is stronger than the control (15) in terms of 
the quantity MA"|. On the other hand, at variance with the estimate (15) above, 
whose right hand side grows exponentially fast in t, the bounds (2.19)-(2.20) in 
[5] involves a right hand side growing much faster in t — specifically, of order 
exp(ci exp(c 2 |t|)) for some constants ci,C 2 > 0. Another difference between [5] and 
our work is that the inequality (15) does not postulate any dependence of e in N 
— on the contrary, e and N are independent throughout the present paper. 

In [14], for each sequence e = e(N) —> 0 as N —> 00 and each monokinetic solution 
of the Vlasov equation (5) i.e. of the form f(t,x, £) = p(t,x)S(i; — u(t,x )) 

Theorem 1.1 gives an asymptotic approximation rate for the convergence of 
the Wigner transform at scale e(N) of p\ N to / in the sense of distributions for 
t € [0,T]. A priori, the time T and the convergence rate depend on the Vlasov 
solution / and on the sequence e(N). On the diagram of section 1.3, this result 
corresponds to the left vertical and bottom horizontal arrow along distinguished 
sequences (e(N), N) over time intervals which may depend on the dependence of e 
in terms of N. Another approach of the same problem can be found in [25]: it is 
proved that each term in the semiclassical expansion as e —► 0 of the quantum IV- 
body problem converges as N -A 00 to the corresponding term in the semiclassical 
expansion of Hartree’s equation. 

On the contrary, Theorem 2.4 provides a quantitative estimate of the distance 
between the solution of the Hartree equation and the first marginal of the solution 
of the quantum V-body problem, for a rather general class of initial data. This 
estimate implies that the mean field limit, i.e. the top horizontal arrow, is uniform 
as e —> 0, over arbitrary long time intervals. This estimate is the quantum analogue 
of the Dobrushin estimate [7] for the classical mean field limit — see section 3. 

The new ideas used in the proof of Theorem 2.4 are 

(a) the use of the quantity M A|, which behaves well with the Toplitz quantiza¬ 
tion, and can be conveniently compared with the Monge-Kantorovich distance with 
exponent 2 on symbols, to which it is obviously analogous; 

(b) an “Eulerian” version of Dobrushin’s estimate, which avoids the traditional 
presentation in terms of particle trajectories as in Dobrushin’s original work [7], 
and can therefore be easily adapted to the quantum dynamics; 

(c) and the adaptation of Dobrushin’s estimate to the V-particle Liouville equation, 
thereby avoiding the need of any quantum analogue of the classical notion of N- 
particle empirical measure. 
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The Eulerian version of Dobrushin’s inequality (item (b) on the list above) sig¬ 
nificantly simplifies the original argument, and allows extending Dobrushin’s in¬ 
equality to Monge-Kantorovicli distances with arbitrary finite exponents (see [21] 
for the original argument for the Monge-Kantorovich distance with exponent 2). 

Estimating directly the Monge-Kantorovich distance between the first marginal 
of the IV-particle distribution function and the solution of the mean field equation 
(item (c) on the list above) avoids using the fact that the IV-particle empirical 
measure is an exact solution of the mean field equation, an important feature in 
Dobrushin’s original approach [7]. This feature is very peculiar to the mean field 
limit in classical Hamiltonian mechanics, and we do not know of any quantum 
analogue of the notion of IV-particle empirical measure which would exactly satisfy 
the mean field quantum dynamics. In other words, the mean field limit in quantum 
mechanics cannot be reduced to the continuous dependence of solutions of the 
quantum mean field equation in terms of their initial data, in some appropriate 
weak topology. 

The outline of the paper is as follows: in the next section, we present items 
(b)-(c) above on the mean field limit for the classical Liouville equation, leading to 
the Vlasov equation. The resulting estimate in Theorem 3.1 below improves earlier 
quantitative bounds of the same type obtained in [15, 22]. The proof of the uniform 
in e estimate in Theorem 2.4 for the quantum mean field limit occupies section 5. 
The properties of the quantity MA"J used in this estimate, stated in Theorem 2.3, 
are proved in section 4. The material on Toplitz quantization, Wigner and Husimi 
functions used in the proof of Theorem 2.3 is recalled in Appendix B. 

The Toplitz quantization is important in the analysis presented here, since it 
connects our main result (Theorem 2.4) with the Dobrushin proof of the mean field 
limit for the classical IV-body problem. For this reason, we dedicate our work to 
the memory of our friend and teacher Louis Boutet de Monvel, in recognition of his 
great contributions to the theory of Toplitz operators. 

3. The Mean Field Limit in Classical Mechanics 

As a warm-up, we first discuss the mean field limit for the IV-body problem 
in classical mechanics, leading to the Vlasov equation. The approach proposed in 
[24, 6, 7] is based on the fact that the phase-space empirical measure of a IV-particle 
system governed by the Newton equations of classical mechanics is a weak solution 
of the Vlasov equation (5). The estimate of the distance between the IV-particle 
and the mean field dynamics obtained by Dobrushin [7] can be formulated in terms 
of propagation of chaos for the sequence of marginals of the IV-particle distribution, 
as explained in [15, 22]. 

The approach proposed below bears directly on the IV-particle distribution, i.e. 
the solution of the Liouville equation (4), and avoids any reference to the IV-particle 
empirical measure. Besides, the core of our argument also avoids using particle tra¬ 
jectories and is based on a computation formulated exclusively in terms of Eulerian 
coordinates. For that reason, this approach can be adapted to the quantum prob¬ 
lem, at variance with the Dobrushin procedure [7], also used in [15, 22]. 

For q, v S T^R^), we denote by n(/z, v) the set of couplings of q et u, i.e. the 
set of Borel probability measures 7r on R d x R d with first and second marginals 

(18) 


7Ti = /.t and 7T2 = v . 
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In other words, 


(19) 


<t>(x)n(dx) + / ip(y)is(dy) 


/ R d 


+ %p{y))n(dxdy) = 

R d xR d J R d 

for each G Cb(R d ). The identity (19) can be used as a definition of the first 
and second marginals of 7 r in (18). Finally, we recall the definition of the Monge- 
Kantorovich distance of exponent p > 1 on V p (Rr): 

_ i /p 

x — y\ p ir(dxdy) 


(20) dist M K, P (/b v ) '■= inf 

7 TGII (/It,!') 


R d xR d 


Theorem 3.1. Assume thatV G Cf(R d ) satisfies (1). Let f m G V p (R d xH d ) with 
p > 1. Let Fn be the solution of the Cauchy problem (4) for the N-body Liouville 
equation with initial data 

jpin ._ 

and let f be the solution of the Cauchy problem for the Vlasov equation (5) with 
initial data f m . For each integer n such that 1 < n < N, let 

F£(t) := J F N (t, dy n+ i ... dy N dp n+1 ... dp N ) G V p { R d ” x R dn ) 

be the n-th marginal of Fn (i.e. the marginal corresponding to the phase space 
distribution of the n first particles). Then 

[p/2] + 1 e A *>‘ - 1 


-distMK.p (f(t)® n ,F£(t))P < 2 p K p \\S7V\\ 


L °° A rmin (p/ 2 ’ 1 ) An 


where K. p := max(l,p — 1) and A p := 2K p (l + 2 P 1 Lip(VH) p ). 

Observe that one obtains an estimate of distMK,p(/(!) 0r \ Ff)(t)) of order 0(1V -1 / 2 ) 
if 1 < p < 2, corresponding to the rate predicted by the central limit theorem. 

Proof. The proof of Theorem 3.1 occupies the remaining part of the present section. 
3.1. The dynamics of couplings. Let 7G , (y™)®-^) satisfy 

(21) T a #n™ = 7 r™ , for each a G & N , 

where 

T,j(xi,tA,... ,XN,fN,Vl,Pl, ■ ■ ■ ,yN,PN ) 

(*£(t( 1 ) 5 £<r(l) ) • ■ * i X<r(N) 5 f,cr(N) 5 Va( 1 ) j V( 7 (l)i * * • ? Va{N) 5 V&(N) ) • 

It will be convenient to use the following notation 

X N ■= (xi,..., Xn) , Ejv := (£ 1 ,..., £n) , 

Yn ■= ( 2 /i, • • ■ , 2 dv), H n := ( 771 ,..., pn) ■ 

Let / be the solution of the Cauchy problem (5) with initial data /*", and let 

N 


( 22 ) 


H 


SW(X N ,E N ) 2 + ^[/]fe)) 

i =1 

be the mean field Hamiltonian. On the other hand, let 

N , N 


U N {Y N ,H N ) := ^2 ll^l 2 + 2 tv S v {Vk — yi) 


k =1 


k,l—l 
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be the microscopic (7V-particle) Hamiltonian. Finally, we denote by {-, -}jv the 
Poisson bracket on (R d x R^)^ defined by 


JV 

3=1 

The following observation is the key to the Eulerian formulation of the Dobrushin 
type estimates, which we shall adapt to the quantum case. 

Lemma 3.2. Let t H > 7T/v(t) G V' 2 ((R d x R d ) N x (R d x R d ) N ) be the solution of 
the Cauchy problem 


dt^N + v) + 'Hn(Yn, Hn),ttn}2n = 0, 

TwLr, = k 1 n ■ 


Then 


7rjv(i) G n(/(t) 0iv ; i'V(t)), for each t G R, 


and 


T a TTjsr(t) = 7Tjv(f) /of each t G R and eac/i a G &n ■ 

Proof. Let f> = Sjv) and = iP(Yn,Hn) G C£°((R d x R d ) N ) be two test 

functions. Then 


d_ 

dt 


d_ 

dt 


JJ (p(X N ,E N )n N , 1 (t,dX N d'E. N ) 


<t>{X n, Sjv)7nv(i, dX N dE N dY N dH N ) 


JJ JJ 4>}n(Xn, 2Ar)7Tiv(t, dX N dE N dY N dH N ) 



pint)} 

N 


<I>}n(Xn, Sjv)7Tjv, 1. (t, dX N dz.N) 


since 

{R-n(Yn , H N ),(j>{X N , E n )} 2 n = 0 . 

The penultimate chain of equalities shows that the first marginal 7T/v,i of ttn cor¬ 
responding to the phase space variables (Xjv,£jv) is a solution of the equation 

d t n N ,i + {H^ /(t)1 ,7rAr a }iv = 0. 

On the other hand, an elementary computation shows that the solution / of (5) 
satisfies 

3 t /^ + {Hf t)] ,/®% = 0. 

Since 71^,1 and f® N are solutions of the same Liouville equation and 

7TA > l(0) = (/ m )® Ar = /(0)® Ar , 

we conclude from the uniqueness of the solution of the Cauchy problem for a trans¬ 
port equation with Lipschitz continuous coefficients that 7Tjv,i(i) = f(t)® N for all 
t > 0. 
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Similarly 


d_ 

dt 


d_ 

dt 


ip(Y N , H N )TT N}2 (t, dY N dH N ) 


4>{Yn , H N )-ir N (t, dX N d^NdY N dH N ) 


{Hn, 4>}n(Y n , H N )ir N (t, dX N dE N dY N dH N ) 


{Hn, iP}n{Yn , H N )wN,2(t, dY N dH N ), 


since 

{H^ /(t)1 (X N , ~ N ), *I>{Y n , H N )} 2N = 0. 

This shows that the second marginal 7 Tjv,2 , corresponding to the phase space vari¬ 
ables (Yn,Hn), is a solution to the same Liouville equation (4) as Fn- Since 

ttn, 2 ( 0 ) = (f in )® N = Fn(0) , 

we conclude that ttn 2 {t) = Fjv(i) for each t > 0, by uniqueness of the solution of 
(4) - 

Finally the time-dependent Hamiltonian 

h / : {X N , Sjv, Y n , H n ) i-a {X n , E n ) +'H N (Y Nl H N ) 


satisfies 


Hence 


Since (ttn 


hf o T a = hf for all cr G S at . 

d t (n N - T a #n N )(t) + {h f: ( tt n - T a #n N ){t)} 2N = 0 . 
F (T #7rAr)(0) by (21), we conclude that 


7Tjv(f) = T a #TTN(t) 


for allt > 0, by uniqueness of the solution of the Cauchy problem for the Liouville 
equation with Hamiltonian h/. □ 


3.2. The Eulerian variant of the Dobrushin estimate. Set 


D p N (t) ■= j - Vj\ p + 16 - Vj\ P Mt,dX N dE N dY N dH N ). 

3 = 1 

We recall that p > 1. Then 

: j |Hf (t)1 (Xv ) S7v),^(|^- 2/j |^ + |e^%n| dn N (t) 
+ j + dn N (t). 


dP p N 

dt 


1=1 


2 N 
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u p y w \x N ,E N ), - - vj\ p +16 - m\ p ) 


I N 1 N 

= M EU&f. \ X j - %| P }2iV + - E{U[/]fe). 16 - %| P }2iV 


= Ti E6 • (*j ^ yj)\xj -yA p 2 - tt E v U[/]6v) • (6 - %)I6 - 2 , 


H N (Y N ,H N ), - - vA p + 16 - Vj r) 


TV N 

= T7 \ X J - %| P }2iV + T7J E - WO, 16 - %r}2iV 


= afE^' ' (y* _a; j)lyj - 


-2 P 


E vn%-^)-(%- 6 )i%- 6 r 2 - 


Therefore 


^ E /(6— Pj) ■ fa - vi)\ x 3 - %r 2 ^ 

j=l 

V N f ( 1 N \ 

-jj E / (6 - vj) ■ vv P [ f] {xj) - — E w (% - y*) 16 - '/.yf ; ' 2 ^.v • 

j=i J V *:=i / 

At this point we use Young’s inequality in the form 

pa^" 1 < a p + (p - l)6 p < max(l,p - l)(a p + 6 P ) 
for each a, b > 0 and each p > 1. Denoting 


K p := max(l,p — 1), 


one has 


dt ~ N 


N r r N .. 

E J (16 - Vj\ p + l*j - yj| p )d7rjv + -j*- E J 16 - %r d7r Jv 


+ E / vv p[/i (*j) - ^ E w (yj - yfc) d7r ^ 

j=l J k =1 

AT N p 

— ‘ 2 KpD p N + E/ vv p[/i ( x j) - 4 E w (yj ~ yfc) d7r ^ • 
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Let us decompose this last term as follows 
1 

N 


1 N , N 

W p[/]Oj) ~ -Vk) = VV p[f] { Xj ) - — J2^V{xj -Xk) 


k= 1 


fe= 1 


1 ^ 

+ ^ 51 ( W “ x k) - W (yj -Vk)) » 


fc=1 


so that, by convexity of the function 2 4 z p on (0, oo) for p > 1, 


1 N 

W p[/] (*j) -]yE W (W “ »*) 


< 2 P “ 


fe=l 

N 


w p[/i(* j) - 77 H w (*j - a ’ fe ) 


fc=l 


+ 2 P_ 


AT 


AT 


^(Vy(xj - x fc ) - W (yj - y k )) 


k =1 


Then, by the same convexity argument as above 


AT 


y^(VF(xj - Xfe) - W (yj - y k )) 


k =1 
AT 


< ^ E i w (^ - **) - - »)!* 


fc=i 


A/ - 


^Lip(W) p ^ u , , 

< —-— 2^ _ Xk > ~ w ~ yfc)l 


N 


k =1 
N 


fc=i 


IV 


Hence 
dD p 


2 p ~ 1 K p Lip(VH) p 


AT 


dt 


N < 2K p D p n + -—~ %I P + l x fe ~ yfe| p )rf7Tiv- 


j,fe=i • 


2 p ~ 1 K p 

aT - 


AT J N 

■ / vv cw(**) (*j -**) 

i=l J k-l 


dirjsf. 


Since 


1 

V 2 


A 2 a 

X! / (I^J “ %l p + kfc - yfc| p )d7Tjv = y] / 1*1 — 2/i| p d7TA < 2D p n , 

;,fc=i ' v ;=i J 


j,k= 1 

the inequality above can be recast as 


^<A p D^ + ^t 


dt 


Wp[/](*i) , r w( x j x k) 


A 


TV 


fc=i 


with 


j=i 

A p := 2AT p (l + 2P" 1 Lip(VH) P ). 


d7Tjv 
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3.3. Controlling the consistency error. Let us examine the last term on the 
right hand side of this inequality. Since 7rjv(t) £ , FN(t)) by Lemma 3.2 


N 


^Vp[f]( x j) jvEWfo Xk ) 


fc=1 


G?7TjV 


JV 


w p[/](*j) - ^ H W (*J - Sfc) 


fc=l 




Observe that this term involves only the factorized distribution f® N . This is the 
“consistency” error in the sense of numerical analysis. In other words, it measures 
by how much f® N fails to be an exact solution of the IV-body Liouville equation 
(4). Eventually we arrive at the inequality 


(23) 


N 


c N 

-£ 

3 =1 ' 


N 


Wp{f}i x j) jy yw( x j x k ) 


fc=i 


f j n p 

< A D P 

dt - P N 


p[f]® N (t,dX N ). 


The last term on the right hand side of (23) is mastered by the following inequal¬ 
ity 1 . 

Lemma 3.3. Let F be a bounded vector field on R d , and p be a probability density 
on R d . For each p > 0 and each j = 1, ..., N, one has 

I P N 


F * p(xj) ~ jtYI F ( x i ~ Xk "> 


N 


N 


k =1 


II p( x m)dx m < jP/ in ( p/2 p)( 2 \\ F \\L°°) P ■ 


The proof of Lemma 3.3 is deferred until Appendix A. Inserting the bound in 
Lemma 3.3 in the right hand side of the differential inequality (23) leads to 


dm 


— < A D p + 2 2p K + ^ 


p N min(p/2,l) ■ 


Gronwall’s inequality implies that 


(24) 


DW) < DU Ole*-' + 2^K r \\VVr L . ^±F e -X 1 . 


3.4. End of the proof of Theorem 3.1. By Lemma 3.2, 

T a #TT N (t) = W N (t) 

for each t > 0 and each permutation er £ ©at. Thus, for each k = 1,..., N, one has 


J(\xk - Vk\ p + |£fc — rik\ p )n N {t, dX N dE N dY N dH N ) 
= J (|ah - yi\ p + 16 - »7i| p )7rjv(i, dX N dZ N dY N dH N ), 


^We thank one of the referees who suggested using this argument instead of the error estimates 
in terms of Monge-Kantovich distances for the law of large numbers obtained by Fournier and 
Guillin [11]. 
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so that 


D p N {t) = J (\xj - yj\ p + |- r]j\ p )'K N {t,dX N d'Z N dY N dH N ) 

r. n 

■ J - Vj \ P + I - T]j\ p )n N (t,dX N dB N dY N dH N ) 


1 

n . 


Denoting 

(25) X% := (xk, ■ ■ ■ ,x N ), E% := (£k,...,£ N ) 

for each k = 1,... N, we set 

tt" (t, dX n dZ n dY n dH n ) := j 7 r N (t, dX™dZ n N dY™dH™). 

By Fubini’s theorem, for each (j) G Ck((R 2d ) n ), one has 

J* 0(-^n? n dX n d^ n dY n dHri) — J* (j>(^X n ^ i^ n )7Tjy(t, dX n dz^N dYjy dH 

= J cj>(X n ,Z n )f(t)® N (dX N dZ N ) 

= j <j>(X n ,Z n )f(t)® n (dX n d~ n ), 


while 


H n ^7Tjy(t) dX n d^ndY n dHn) — (f){Y n , R n )7TjV (t, dX jy d^j^dYj^fdH jV ) 

= J <KY n ,H n )F N (t,dY N dH N ) 


so that 
Therefore 

(26) 


= J <j>(Y n ,H n )F%(t,dY n dH n ), 

7r“ (£) G n(/(t)®", F"(f)) for each t> 0 . 

1 /' " 

W) = - / Ed^' ^ ^ + l& - vA P )*N(t,dX n dZ n dY n dH n ) 
n J ,_i 


i=i 


> -distMK,p(/(*)® n ,i^(t)) P . 
Thus, inequality (24) implies that 


(27) 


i dist M K, P (/(t)®",i^(i)) p <D^(0)e Ap 


Ar,t 


■2 2p K p \\XV\\ p L 


[p/2] + 1 e Apt - 1 


L °° ATmin(p/2,l) A 


for each t > 0, each AT > n > 1 and each tt 1 x G n , F‘™) satisfying 

T a for all a G 6jv . 

Finally, choose the initial coupling of the form 

^ := (f n f N (dX N dZ N )S {XN ^ N) (Y N ,H N ), 
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i.e. 

n N := D#(/ ira )® JV , where D : (X N , ~ N ) ^ (X N , Z N ,X N , Z N ). 

Then 

1 N r 

(28) D p n ( 0) = - ^ / (\x k -y k | p + |& - rj fc |»’)7r}y(<iX Jv da Jv dYi v dffjv) = 0, 

iV j=i ^ 

since 7r“ is supported in the diagonal 

{(X/v, Z N , Y n , H n ) s.t. X N = Y n et Sjv = Hn} . 

Inserting this last piece of information in (27) leads to the inequality in the state¬ 
ment of Theorem 3.1. □ 

Remark 3.4. In the argument above, we have not checked that f(t)® n and Fff(t) 
belong to V p ((H d ) n ) for all 1 < n < N and all t > 0. This verification is needed on 
principle, because the Monge-Kantorovich distance of exponent p is defined on the 
set of probability measures with finite moments of order p. For each p > 2, set 

M p(t) ■= Jj (W P + |£| p )/(t, dxdf ), t > 0 , 

where f is the solution of (5). An easy argument based on Young’s inequality and 
the mean value inequality for XV shows that 

M p (t) < Aip(0)e (p - 1)(1+2 Lip(w))t , t > 0 . 

In particular f(t ) £ ^(R^ x R d ) for all t > 0 provided that f ln £ ^(R^ x R d ). 
Therefore /(f)®” £ V p ((R d x H d ) n ) for each n> 1, and the inequality in Theorem 
3.1 implies that F(t )" £ ^((R^ x R d )) for all t, > 0 and each n = 1,..., N, since 
distMK,p(/(i) 0 ",^(OAr) < 00 • 


4. Proof of Theorem 2.3: Properties of MI <| 

4.1. The general lower bound (14). The lower bound (14) can be viewed as a 
variant of the uncertainty principle. More precisely 

Q*Q + P*P = (Q + iP)*(Q + iP) + i{P*Q - Q*P ) > i{P*Q - Q*P ). 

On the other hand, 

(■ P*Q ^ Q*P) = ~ie (div Xl (xi) + div X2 (a:2)) = - 2 idel^ 2 , 


so that 

Q*Q + P*P > 2del?, 2 . 
Therefore, for each R £ V(S) 2), one has 


trace r -, 2 ((Q*Q + P*P)R) = trace,-,, (R 1/2 (Q*Q + P*P)R 1/2 ) 

> 2de trace*,, (R) = 2 de . 
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4.2. Upper bound for Toplitz operators. The goal of this section is to prove 
statement (1) in Theorem 2.3. 

In the course of the proof, we shall need the following intermediate result. Ev¬ 
erywhere in this section, the identification R 2d B (q,p) >-> q + ip G C d is implicitly 
assumed. 

Lemma 4.1 (Construction of quantum couplings). Let e > 0 and let p\ and p\ 
be Toplitz operators at scale e on L 2 { R d ) with symbols (2Tte) d p\ and (2Tre) d p 2 , 
where pi,p 2 G ^(R 2 ^). For each coupling p G II(/xi,p, 2 ), the Toplitz operator 
OPf ((27re) 2d /x) belongs to Q(p\,p 2 )- 

Proof. Set p e j = OPf ((27 xe) d pj) for j = 1, 2. We recall that f) = L 2 (R d ). 

Since p is a, probability measure on C 2d , the operator OPf ((27 re) 2d p) belongs to 
T>(Sj 0 S)) by (47). Next, one has 

O ((2ire) 2d p) = [ \(zi,z 2 ),e){(zi,z 2 ),e\p(dzidz 2 ), 

J c 2d 

where the wave function \(zi, z 2 ),e) is defined in formula (45) in Appendix B. In 
this section, we use the bra-ket notation, also recalled in Appendix B. Observe that 

\(zi, z 2 ), e)((zi, z 2 ), e| = \zi,e)(zi,e\ 0 \z 2 ,e)(z 2 ,e\ , 

so that 

tracer#((A 0 7#)|(zi, z 2 ), e)((zi, z 2 ), e|) 

= trace# (A| z lt e) (zi, e|) trace#( \z 2 , e)(z 2 ,e\) 

= trace#(A| z lt e) (zi, e\){z 2 , e\z 2 , e) 

= trace# (A|zi, e){z\,e \), 

for each A G C(S)). Therefore, 

trace#®#((A 0 7#) OPf ((fh:e) 2d p)) 



trace#®#((A 0 7#)|(,zi, z 2 ), e){{z\, z 2 ), e\)p{dz\dz 2 ) 



tva.ce^{A\zi,e){zi,e\)p(dz\dz 2 ) 


= / trace#(A| 2 i, e){zi 1 e\)p\(dz\) 

Jc d 

= trace#(AOPf (( 2ire) d pi )). 

Likewise, for each A £ £(jj), 

trace#®#((7# 0 A) OPf ((27re) 2£ V)) = trace#(AOPf ((2ne) d p 2 )) = trace#(Ap 2 ) ■ 
The conclusion immediately follows. □ 

Proof of Theorem 2.3 (1). For each p G II(/zi, p 2 ), by Lemma 4.1 and the definition 
of Miq , one has 

MK^pi, p 2 ) 2 < trace#®#((Q* Q + P*P ) OPf ((2tt e) 2d p )), 


where we recall that 


( 29 ) 


Qlf(x!,X 2 ) = (xi - x 2 )lp(xi,x 2 ) , 
Pfj(xi,x 2 ) = -ie(S7 Xl - \7 X2 )ip(xi,x 2 ), 
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so that 


Q*Q 4 >{Xi,X 2 ) = \xi - X2\ 2 1p(xi,X 2 ) , 


and 

P*Pip(xi,x 2 ) = -e 2 (div Xl -div X2 )(V Xl - V X2 )ip(xi,x 2 ). 
In other words 


Q*Q + P*P = OPf (\ qi - q 2 1 2 + \ Pl - p 2 1 2 ) - 2delsj 2 

— see formula (48) in Appendix B, with f{yi,y 2 ) = \yi — y 2 \ 2 , and A/ 
formula (55) in Appendix B, one has 

trac e mSi ((Q*Q + P*P) OP^((2^e) 2 V)) 


(30) 



( 1*1 - x 2 \ 2 + |Ci 


f 2 \ 2 ) y(dxidfidx 2 df 2 ) + 2 de. 


4d. By 


Therefore, for each y G H(yi,y 2 ), one has 

MKl(pi,p 2 ) 2 < [ (\x! - x 2 \ 2 + |£i - £, 2 \ 2 )pi{dxid£ >1 dx 2 d£, 2 ) + 2de. 

J( R d ) 4 


Observing that the left hand side of this inequality is independent of y G II(/Lii, ^ 2 ), 
we conclude that 

MK^pi,p 2 ) 2 < inf [ {\xx-x 2 \ 2 + \£,i-^ 2 \ 2 )y{dxid^dx 2 d^ 2 )+2de 

•n-en(/ii,p 2 ) J(rj)4 

= distMK, 2 (Mi, M 2) 2 + 2 de 


which is the sought inequality. 


□ 


4.3. Asymptotic lower bound for MK J. The core of the argument leading 
to the lower bound in Theorem 2.3 (2) combines Kantorovich duality with the 
convergence of Husimi functions to Wigner measures. 

Proof of Theorem 2.3 (2). Let a, b G C^R^) satisfy 

(31) a(aq,£i) + b(x 2 ,&) < \x, - x 2 \ 2 + |£i - £ 2 | 2 for all x 1 ,x 2 ,£i ,&G R d . 
Hence 2 

(OPf (a) <g> If> + If, ® OPj (6)) = OPf (a <E> 1 + 1 <g> 6) 

< OPf (c) = Q*Q + P*P + 2delsj 2 ■ 

Thus, for each R e G Q{pl,p 2 ), one has 

trac e m s,((Q*Q + P*P)R e ) 

> trace i 5 ®j 5 ((OPj’(a) <8)1?,+ Is, ® OPf ( b))R e ) - 2 de 

= trace^,(OP^(a)pi) + tracer(OPf (b)p e 2 ) — 2 de . 

Taking the inf of the left hand side as R e runs through Q{p \, p 2 ), one arrives at the 
inequality 

(32) MK\{p\,p’f) 2 >tvace?,(OPj (a) p{) + tra,ce?,(OPj (b)pl) — 2de. 


2 Denoting a(g)l, 1 <g)fo and c the functions ( x 1 +1, X2+2) >—> a(x, , £1), (xi, £1, X2+2) H> 6(^2,£2) 
and (xi,£i, X2, £2) >-A |xi — 121 2 + |£i — C2I 2 respectively. 
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Next the traces on the right hand side are expressed in terms of the Husimi functions 
of p\ and p| by formula (54): 

tracer (OPf(a)pi) = JJ a(x 1 ,^ 1 )W e [p e 1 ](x 1 ,^ 1 )dx 1 d ^ 1 , 


tracer (OPf (6)p|) 


b{X2,&)W £ [p2](X2,ti2)dX2dt;2 ■ 


Hence 

MKl{p\,p\) 2 > JJ a(x 1 ,^ 1 )W e [p e 1 ](x 1 ,^ 1 )dx 1 d^ 1 

+ JJ b(x 2 ,&)W e [pl](x 2 ,&)dx 2 d£ 2 - 2 de. 

Taking the sup of both sides of this equality over all a, b G Cb( R d x R d ) satisfying 
(31) shows that 

MKfolptf 


> sup 

a<g)l+l<g)b<c 

a,beC b (R d XR d ) 


aW e [p{]dxid£i 


bW e [pl]dx 2 d& 


— 2 de 


= dist M K,2(W / e[Pi], W e [p e 2 ]) 2 - 2 de, 


where the last equality follows from Kantorovich duality (Theorem 1 in chapter 1 
of [31]). This gives the first inequality in Theorem 2.3 (2). 

Since the sequence of Wigner transforms of the density matrices p e j satisfies 

W e [pj] Wj in <S'(R d x R d ) as e — > 0 

for j = 1,2, one has 

W f [pj] Wj weakly in the sense of measures on R ci x R d 

for j = 1,2, by Theorem III. 1 (1) in [20]. By the first inequality in Theorem 2.3 
(2) already established above and Remark 6.12 in [32], 

lim MK%(pi , p %) 2 > limdist M K, 2 (W e [pi], W e [p e 2 }) 2 > dist M K,2(ub, w 2 ) 2 ■ 

e—>-0 e—>-0 

This concludes the proof of Theorem 2.3. □ 


5. Proof of Theorem 2.4 
The quantum 7V-body Hamiltonian is 

N 1 

We,N ■= ^ — 2 e A 7 + ’ 

i=l 

an unbounded self-adjoint operator on -fjjv := $)® N = L 2 ((H d ) N ). We denote 
by A j the Laplacian acting on the variable Xj, and by Vp- the multiplication by 
V(xj — Xk )■ The IV-body von Neumann equation for the density matrix is 

(33) ied t p e ,N ■= [Tde,N,Pe, n] , Pe,w| t=0 = P^N ) 

with initial data G 2?Cfjjv). 

On the other hand, for each p G T>{9 j), we consider the mean field Hamiltonian 

He :=-±e 2 A + K p , 
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where V p designates the operator defined on S) as the multiplication by the function 

x i-A J V(x — x')p(x', x')dx' =: V p (x). 

The Hartree equation for the density matrix p e is 

(34) iedtpe = [Hf ,p e \, Pe\ t=0 = pT ■ 

By a straightforward computation, the solution p e of the Hartree equation (34) 
satisfies 


(35) 

where 


ied tP f N = \H p ' N ,pf N ] 


®N _ / in\®N 
J e t =0 ^ > 


Kn ■= Hf 0 I 


f)N -1 


■ I ® Hf ® L 


f)N -2 




f)N -1 


Hf 


5.1. The quantum dynamics of couplings. Let £ Q((pl n )® N , P™n)-> and 
define t i-a R e ^(t) G a) to be the solution of the Cauchy problem 

(36) ied t R e ,N = [H^V <g> I^ N + I^ N ® W e ,N, Rc,n] , Re,N\ t=0 = R™n ■ 

For each <r G &n, we denote by T a the unitary operator on ^ 2 a defined by 

5 • • * 5 Xn ; Vli • • • ) 2/a) 4 /(x ct -i 5. * • * j Xfj—i (A) 5 2/cr — qi) 5 • • • ) 

for each 4/ G Jo 2 a- We shall henceforth assume that 

(37) %R™ n T* = R™ n , for each a G ©a- 

Lemma 5.1. For eac/i t > 0, one has R e ,A (t) G Q((p e (t))® N , p e . Ar(t)). Moreover 

(38) T a R t ^N(t)Ta = i? £i Ar(t) for each a G ©a and t > 0 . 


Proof. By definition 

Re,N (t) = Ue t N{t/e)Rl n N U et N(t/e)* , 

where 

-jjU et N(t) = ® + /r, N ® 'H €> A)b r e ,A(t), C’ £ ,at( 0) = /f, 2JV . 

Since V ® © TF.a is self-adjoint, t/ e ,A(t) is unitary for each f > 0. 

Therefore 


Re,N(t) = R e ,N(t)* > 0 and tracer, 2N (R e ,N(t)) = trace ,-, 2N {R™n) = 1 


for each t > 0. 

The marginals of the density 7 ? £i at in the product f) 2 A = f)A ® f)A are defined 
by analogy with (18)-(19) in Definition 2.1: 


(39) 


Re,N,i G T>{?)n) and tracer , N (AR £j a,i) = trace,-, 2JV ((A ® I^ N )R e ,N), 
Re,N ,2 G V(Sj N ) and trace,-, N (AR £j a, 2 ) = tracer, 2N ((Isj N ® A)R e , N ), 


for each H G C(Sjn)- 

Let A G C(S"jn) be such that 




l=i 


G C($j]\r )> 
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then 

ie<9 t tracer(AR^ jv.i) = *e<9 t trace$ 2JV ((A 0 Isj N )R e , N ) 

= - trace^ 2JV ([H^' JV ® I^ N + J fljv 0 H e ,N, (A ® -T^jv)]^ e,jv) 

= — trace^ 2N (([H^jy, A] ® Isj N )R e ,N) 

= ~ tracesj N ([U^ N , A]R e<N>1 ) 

= trace ftjv (A[H^ e jV , tf e ,iv,i]) 

and 

p I _ p in _ / in\®N 

Re,N,l\ t=0 - ,JV ,1 ~ (P ) 

By uniqueness of the solution of the Cauchy problem (35), one concludes that 
Re,N,i{t) = p e (t)® N , for each t> 0 . 

Similarly 

ied t tracer (AR etNt2 ) = ie9ttrace ft2jv ((/ i5jv 0 A)R e , N ) 

= - trace ft2JV ([H(’ e JV ® I^ N + I^ N 0 U eyN , {I^ n ® 2 l)]i? e>iV ) 

trace^ 2 iv ((-^)ojv O [Hey, -tl])_/? e yY) 

t race^ ^ ([% €j a ? A\ Re,N,2) 
trace^ N a? Re, a, 2 ]) 

and 

p I _ p in _ in 

-fte,A,2| t=0 — -tt e ,A,2 ~ Pe,N ■ 

By uniqueness of the solution of the Cauchy problem for the von Neumann equation 
(33), one concludes that 

Re, a, 2 (t) = Pe,N(t), for each t > 0 . 


Finally, let cr £ &n', 

iedt(TaRe,N(t)T a *) = Ta[H-£* N 0 Is) N + If) N 0 He,A, -Re,A]7^ 

= [7^-(Hg'jy 0 + Ifi N ®'He,N)'Tai'T a Re,N'Tv] 

= ® / fiK +4^ 0 T a -u ttN T:,%R e , N T;] 

= [H^ 0 + Isi N 0 W e ,A, TaRe,NT*} 

where r CT is the unitary operator on 1 )jv defined in (9), since one has obviously 
r n H^r; = and T ff H £ , N r; = ?4 ,a - 

On the other hand 


%Re, N (o)r; = r a Ri n N (o)T; = r? n = Re, N ( o), 


so that, by uniqueness of the solution of the Cauchy problem (36), 

T a R e ,N(t)T* = Re, N {t) for each t > 0 and each a e &n ■ 


□ 
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5.2. The Dobrushin type estimate. Set 


N 


D e , N (t) := trace — y~]( Q*Qj + P* Pj)R e , N (t) . 


j=i 


We compute 


dD , 


dt 


N 


e,N =-trace + Pj Pj)[^N ® If>N + If>N ® Pe,N, Re,N] 


3 =1 


= — trace | i 
e 


= — trace | i 
e 


H— trace | i 
e 


N 


® ® ^e,AA -jy YMQj + Pj Pi) 


3 =1 
AT 


jy 0 /qJV — He,iV 0 /fjjv , ^E(^+a*a) 


i=i 

N 


P-e.N ® Aw + 4 n ® Re.Ni Y^QjQi + PjPj) 


i=i 






Re 


,N 


,N 


N 


Consider the second term on the right hand side of this equality. One has 
[He,JV ® A W + A W ® H e ,JV, (Q*Oi + P*Pj)\ 

= [~h e (A ® A + A ® A)i QjQj]j 


2 " 
N 


+ ^ EllUfc ® A + 1$ ® A'fe’ A Rj]j 


fc=i 


where the index j on the brackets in the right hand side indicate that the corre¬ 
sponding operators act on the variables Xj et yj. In other words 

A j = I% U ~ 1) ®A®I% {N - j) 

for each operator A on fy. 

Obviously 

[-A 2 (A ® A + A ® A). QjQj] 

= - A 2 ( d A [Vxj, QjQj]+ [V x ,, Q*Qj] ■ V Xj ) 

- y (2 ''liv, ; [V,,, QjQj] + X yj , QjQj] ■ V v ,) 

= -ie(P*Qj+Q*Pj). 

On the other hand 

[V jk ® Isi + I* ® V jk ,P*Pj] =[V jk ® hi + A ® V jk ,P*]Pj 

+ Pj [Vjk ® hi + A ® Vjk, Pj] 

=ie(yV(xj - x k ) - W(j/j - y k ))Pj 
+ ieP* (W ( Xj - x k ) - W(% - y k ))). 
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Therefore 


- trace | i 
e 


N 


H e ,N <S> Ifi N + If) N 0 He, AT, ^y 'ZiQjQj + PfPj) 


3 =1 

1 


R, 


e,N 


N 


= trace — y/(P* Qj + Q*Pj)R e , N 


i=i 


- trace — ^ P* (W {xj - x k ) - W {yj - y k ))R e , N 


1 <j,k<N 


- trace W (% “ Vk))PjRe,N 


l<j,k<N 


By the Cauchy-Schwarz inequality 

| trace (( P*Q 0 + Q*Pj)R e , N ) | < trace ((P/P, + Q^Q^R^n) , 

while 

trace((P/(VF(xj - - 2 /fc)) + (VV r (a; i - x k )-VV(yj - y k ))Pj)R^ N ) 

< trace(P/PjP ei jv) + trace(|W(a:.,- - x fc ) - - 2 /fe)| 2 P e ,Ar) 

Eventually 


- trace | i 
e 


N 


H e ,iV ® If) N + /$at ® He,AT, ^y ^^(QjQj + PjPj) 


3 =1 


He 


AT 


AT 


< trace — ^(P/P,- + Q*Qj)R^ N + trace — ^ P/P,P e ,jv 


i=i 


iV 2 


l<j,fc<AT 


+ trace |VV r (a; i - a;*,) - W(j/j - yfe)| 2 P e ,iv 


i <3,k<N 


< 2P e ,jv + Lip(VF) 2 trace ^ |(£j - x fc ) - (% - yfc)| 2 P e ,jv 


l<j,k<N 


< 2D €)JV + Lip(W) 2 trace ^ (kj ~ %'| 2 + kfe ~ yfc| 2 )A,JV 


l<j,fc<JV 


= 2P ejA r + Lip(vy) 2 trace ^ ^ |x; - yi\ 2 R e , N j 


< (2 + 4Lip(W) 2 )£> e|JV! 
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so that 


dP e ,N 

dt 


< (2 + 4Lip(W) 2 )D £iJV 


H— trace | i 
e 


N 


H P t *N ® Isj N - Ue,N 


3 =1 


P 


•e,JV 


Next observe that 


N 


H P e e N ® If) N ~ Ue,N <8 If) N , 


i=i 


1 * 1 

= ^ em>.(**)>^ E [ y (^ - *o, (p^+ p?m 


3 =1 


l<fc<Z<iV 
N 




i=i 




l<k<l<N 




iV 2 


1 <k<l<N 


i * 

= ^E([^ + W. (*i)v^-D 


i-i 


iV 2 


E ([v(* fc - xO,p fc *]P fc + Pfc [v(® fc - *0,^]) ■ 


l<fc,Z<N 


Moreover 


(*i), P/]Pi + P*[V Pt fo), Pi] = ie(VV Pe (x,)Pi + P/W„. (Xi)), 


while 

[V(a: fc - *,), P fc 1 Pfc + P fc *[V(x fc - x t ), P k } 
= ie(VV(x k - Xl )P k + P£VV(x k - x t )). 
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Therefore 


8 I^ N — 'He,iv <8 Ifj N , 


= E( w «..feW + p;vD.fej) 


+^2 E (W(* fc -*,)Pib + PfcW(* fc -a:,)) 


1 <k,l<N 


l N ( 1 N \ 

= E m E w (^ - *0 p j 


i N ( i N 

-iv£ p ; w *fo) - af E 


i / i 

- trace I i H^ jV ® H, N - He ,a 8 H..v YliQjQj + P j P j) H € ,a 

= -ij 2 trace ( ( VV IA<Xl) - T7$I VV, (*J -**)') P i^,A 


- T7 X] trace ( p i ( VV P* (*i) ~ 4 X) W (*J _ *fc) ) 


1 V x IV 

< T7 E trace VVp 6 (xj) - — ^ W- x fc ) H £ 


1 N 

— trace (P*PjR e>N ) , 


so that 


1 / 1 1V 

- trace I i H^ N 8 - H e , a 8 , j^EWi+W ^ 

1 a / L a 

< A, A + T7 E traCe VV P. (*j) (*i - Xfc) He,A 
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Therefore, the differential inequality satisfied by -D 6 ,at is recast as 


dD ( 


dt. 


< (3 + 4 Lip(VF) 2 )ZA, 


JV 


JV 


A 


+ ^X trace 

3 = 1 

= (3 + 4Lip(W) 2 )L> e ,jv 

N 


VVp e (Xj) - X W ( X J ~ Xk "> 


k =1 




.A 


1 \ ' 

— > trace 


iV 


i=i 


A 


VF Pe ( Xj )--E VF (xj - x k ) pf N 


k=1 


The last inequality comes from the fact that the multiplication by 

.2 


VF Pe (+,) - ^ X W ~ Xk> > 


A 


N 

fe=l 

acts on the variables variables sq,... , aqv only, and that the first marginal R €) n,i 
de jv is pf ^ by Lemma 5.1. 

Let us finally estimate the term 

i 2 

1 

N 


A 


1 \ " 

— > trace 

N ^ 

3=1 

By a straightforward computation 

trace 


A 


W Pe (:r j )--£vV(x j -x k ) pf N 


k=l 


N 


VF Pe (xj) - - X VFfo - *0 pf N 


k =1 


1 * 

VF(a;.,- — x')p e (t, x ', x')dx' — — VF(aq — aq) 


fe=i 


A 


JJp e (t,aq,aq)dX';v, 


t=i 


Set 


/e(M) ; = Pe(t,X,x) , 

where p e is the solution of the Hartree equation (34). Since p e (t) € T>(Sj), then f e (t) 
is a probability density on R d for each e > 0 and each t > 0. Applying Lemma 3.3 
with p = 2 and p = f e (t, •) shows that (see formula (44) below) 

1 2 


trace 

Therefore 


1 


A 


VF Pe {Xj) - — X W ( X 3 - x k ) 


k =1 


pf w — yllvF||| 0 


^W<AA,aW + ^||VF|||» 


with 


A := 3 + 4Lip(VF) 2 . 
Gronwall’s inequality implies that 
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5.3. From to MA"|(p e (t) 0iv ,/9 e ,jv(t)). Let us denote by R,f N the marginal 
density of P £> at corresponding to the n first particles. In other words, for each 
A € £(f) 2 n), we set 

trace^ 2 „ (AR? N ) = traces) 2 „(S(A ® Isj a( „_ n) )S*R £iN ) 

where 

S&(xi,y 1 ,x 2 ,y 2 , ■ ■ -,x N ,y N ) ='S/(x i,.. .,x N ,yi,.. .,y N ) ■ 

We claim that 

(40) R? >n G Q(pf n ,Pe, N )- 

Indeed, for each B € £(Sj n ), by Lemma 5.1, one has 

tracej 5 2n ((S 0 IsjJR^n) = trac e% 2N (S(B 0 ® Isj 2(N _ n) )S*R etN ) 

= trace^ 2JV ((S 0 /$„_„) 0 I# N )R eiN ) 

= trac esj 2N ((B 0 I^ N _ n )R e ,N,i) 

= tracej5 w ((S 0 )pf iV ) 

= tracej5 n (Spf n ), 

while 

tracer 2n ((/f,„ 0 B)R™ N ) = trace fl2JV (S(/ fln 0 B 0 I^ (N _ n) )S*R t<N ) 

= trace# 2jv ((/ fljv 0(B0 I^ N _ n ))Re,N) 

= tracer, 2iY ((P 0 Isj N _ n )R e ,N, 2 ) 

= tracer, w ((A 0 h-, N _ n )pe,N) 

= trace# n (Bp™ N ). 

Observe further that, for each j = 1,..., TV, one has 
tracer ((Q*Qj + P*P,)P £ ,iv) = trace* 2JV (T CT * (QJQi + P?PJT^n) 

= tracesj 2N ((QiQi + P,* P^T^R^nV.) 

= tracer ((QJQi + P*Pi)P £>J v), 

where a, is the permutation of {1,..., IV} exchanging 1 and j and leaving all the 
other integers invariant, because of the identity (38) in Lemma 5.1. 

Hence, for each IV > n > 1, each e > 0 and each I > 0, one has 

D e>N (t) = tracer,^ ((Q*Qj + P*Pj)R e ,N(t)) for each j = 1, ..., N, 

1 / n 

= - tracer, 2iv I ^{Q*Qj + P*Pj)R^ N (t) 


= ~ tracer, 2n ( 'f^iQjQj + P*Pj)R™ N {t) ] . 


U = ! 


Because of (40), this implies that 


D e , N (t)> -MK* 2 (p e (t)® n ,pZ N (t)). 

n 
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for each N > n > 1, each e > 0 and each t > 0. In particular, 

Qoo pAt _ 1 

(41) MK* 2 (p e (t)® n ,Pl N (t)) 2 <nD e , N (0)e At + -^U W|||oo^— . 

This inequality holds for each R" l N G QM n )® N ,Pt N ) satisfying the symmetry 
condition (37). 


5.4. Upper bound for D £> at( 0). Let Q™ N G Q{(pl n )® N ,P £ (jv)> an d se t 

pm_ ^ \ ' - 7 “ y^)in *7^* 

^e.AT •— 'ffVe.iV'a • 

crG©iv 

By construction, G T>(S)n ® $)n) and satisfies (37). On the other hand, for 
each A G £(S)n), one has 

tracej7 N 0jj JV ((/i5 w ® ^4)i? e ,Af) ^ ' tracer,0 A)TJrQ e ,jvU ) 

crG©iv 

= (/y? ^ ' t ra cej7 W 0f )N (7) T (Ljjjv (8> ^4)7 ^tQ £ jv) 

crG© ;v 

<xG©iV 

= ^T trace^ (r^T^p’^) 

cr^&N 

= ^7 traCe ^JV (^4 T crP £ "jV r CT) 

<re©jv 

= tracer (^p £ » 

since satisfies the symmetry condition (10). By the same token 

trace j5jV 0 ftjv (( £ 4(8) I?, n )R™ n ) = trace fliV (A(p™)® JV ) 

for each A G £(f)jv). Hence 7?“^ G Q({p\ n )® N , p 1 ^)■ 

On the other hand 
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v i ^ +wo r — ^ Bw*+wo 


because 


for each a £ &n- 

Inserting this in (41) shows that 

i o p At _ 1 

-MK%(p e (t)® n ,p™ N (t)) 2 < — ||VF|||oo—-— 


+ — tracer,,, ^{QjQj + P;Pj)Q% 


for all Q l e n N £ Q{(pT)® N , P™ N ) with N >n> 1. Minimizing over Q l ™ N leads to 

I 8 p At — 1 p At 

-MKl{pSr n ,Pl N {t)f < j^llWHioo—^— + W MKI((p?)® N ,piy) 

which is precisely inequality (15). 

Next, in the special case where p l e n = OPj ((2ne) d p,\ n ) with p l e n £ V 2 {C d ), 
while ply = OPj((2ne) dN piy) with p,™ N £ V 2 {(C d ) N ) satisfies (10), one has 
( p m)®iv = OPf((27re) dJV (^ n )® JV ), and we deduce from Theorem 2.3 (2) that 

MK^(p™)® N ,piy) 2 < distMK,2((Mr) 0iV ,^?iv) 2 +2Nde. 

Inserting this in (15) shows that 

i o p At _ i 

-MKi(p e (tr n ,pi N (t)) 2 <^||vy||i„— 


■dist M K,2((Mr)® JV ,^yv) 2 + 2de) 


which is (16). This concludes the proof of Theorem 2.4. 

Remark 5.2. Let p e be the solution of the Cauchy problem for the Hartree equation 
with initial data p l e n £ T>(Sj) such that 

trace((|x| 2 — e 2 A )p l f l ) < oo . 

Then, for each e > 0 and each t > 0, one has 

trace(|a;| 2 p e (f)) < (trace(|a;| 2 p* n ) + fe*(trace(—e 2 Ap“) + 21|V||) < oo , 

as a consequence of the differential inequality 

^ trace (\x\ 2 p e {t.)) = (trace ([|eA, |a;| 2 ]p e (t)) | < trace ((|a;| 2 - e 2 A) p e (t)) 

and of the energy conservation 

^trace(—e 2 A p e (t)) + jj V(x — z)p € (t, z,z)p e (t, x,x)dxdz^ =0. 

In particular, for each t > 0 and each e > 0, one has 

trace((|a;| 2 — e 2 A )p e (t)) < oo . 
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Appendix A. Proof of Lemma 3.3 
Without loss of generality, we treat the case j = 1. Define 
V{z) := F * p(x i) — F(x i — 2 ), 2 : € R d . 

A.l. The case where p > 2 is an even integer. Start with the decomposition 


f 1 A, 

/ F-k p{x 1 ) -—^2,F{x 1 - x k ) 
J V fe=l 


AT 


N 


N^ V{ - Xi) 


Q p(x m )dx 

m =1 
P JV 

p(x m )d 2 


1 

ATP 


E 


.P/2 


m=l 

N 


/ nv(%) ■ ^(^(o) n p{xm)dx m . 

j,fce{i,...,Ar}{ 1 .--.p/ 2 > 1=1 m=i 

For each pair of maps j,k&{ 1,..., define the map {1,..., N jf 1 ’ --^! 

as follows: 


g{l) = j(l) if 1 < l < p/2 , g(l) = k(l - p/2) if p/2 + 1 < Z < p. 

Conversely, given g G {1,..., one reconstructs the maps j,k by the 


formulas 

Hence 


j(0 = ff(0> k{l)=g(l+p/ 2), l<l<p/ 2 . 

N 

p{x m )dx 


A 


1 

NP 


F * p(x 1 ) - — F ( X 1 “ x k) 

k —1 m =1 

* p/2 A 

e /n V{Xg(l)) ' V{,Xg(l+p/'2 )) | | P{x m )dx r) 

g£{ 1,...,A}! 1 .•■■.*>} (=1 m=l 

^ ,. P/ 2 A 

= T7F E / IT V (*9(0) ' V (*fl(i+p/2)) LI P( X m)dXr, 


gts*,” *=i 

/. P/2 


ra=l 

N 


+ j^ e / n v ( :r ffw) • v (^d+p/ 2 )) n 


geM p" i=i 


m=l 


where 

Afj’r := {(/ G {1,..., s.t. #p _1 ({m}) ^ 1 for each m = 2,..., A} , 

Obviously, for each m G one has 


.p/2 


A 


<(2||A|Uoo) p . 


Wy( x g(i)) 'V( x g(i+p/ 2)) P{ x m)dx 7 , 

Z=1 m=l 

Next, for each g G 5^, dehne 

lg ■■= min{Z = 1,... ,p s.t. p(Z) > 1 and #5 _1 ({5(0)}) = !} , 
l g = l g +p/ 2 if l g < p/2 , and = l g — p/2 if l g > p/2 . 


( 42 ) 
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Thus, if g £ S p Nl one has 


.p/2 


N 


• V(. X g(l+p/ 2 )) n p(x m )dx rl 


1=1 


m =1 


\ 


H V(* fl(I) ) • V(a; g ( p /2+z)) 


i l<^<P/2 

\l<£{lg,ig} 


V(fl 




]J[ p(x m )da 


J 


1 <m<N 
m 7 ^9(1 g) 


■ J V i x 9(l 9 ))p( x g(l 9 )) dx 9(l 9 ) = 0 - 


since 


J V( x g(i 9 ))p( x g(i 9 ))dx g (i g ) = J {F*p(x 1 ) - F(xi - z g (i 9 )))p(z g (i 9 ))ete g p 9 ) 
= F*p( Xl )~ J F{xi - x g{ i g) ))p{x g{ i g) )dx g (i g ) =0. 


Hence 


AT 


AT 


np(^ m <^(2iiFiuoo f . 


m=l 


iVP 


(43) f F* p(xj) -jr ^2 F ( x i - x k) 

J v fc=i 

Now 

= /v 2p - #s*, 

and we next compute #Sff. 

For each element g £ 5%-, there are TV — 1 choices for g(Z 9 ) £ {2,7V}, where l g 
is the integer defined in (42). For each such choice, the restriction of the map g to 
{1,... ,p}\{l g } takes its values in {1,..., N}\{g(l g )} and can be chosen arbitrarily 
among the maps from {1,... ,p} \ {l g } to {1,..., TV} \ {g(l 9 )}. Hence 


so that 


Hence 


= (N — 1)(TV — l) p_1 = (TV — 1) 

#M p n = N p — (TV — 1 ) p . 

#M p N 


= 1-11-1)' < Z . 

Np \ TV j ~ TV 


p 


Inserting this inequality in (43) shows that 


(44) 


N 


F * p(x i) - — ^2 F ( X 1 - x k) 


k=l 


N 


m= 1 


II p(x m )dx m < -^(2||F|| L =o) p . 


A.2. The case where p > 2 is not an even integer. Write p/2 as 
p/2 = (1 — 9)\p/ 2] + 9{{p/ 2] + 1), with 0 < 9 < 1, 
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where [z] denotes the largest integer less than or equal to 2 . Using Holder’s in¬ 
equality and (44) shows that 


F * p(x 1 ) - F ( Xl ~ Xk ) 

| ipYl — 1 

< (2[ P /2]) 1 "((2b/2] + 2)^ (2||f||Loo)p 


N 


F[ P(Xm)dx r 


P 


N 


N 


< 


2 [p/ 2] + 2 
N 


m\F\\ L ~r. 


A.3. The case where 0 < p < 2. Using Jensen’s inequality and inequality (44) 
for p = 2 shows that 


/ F*p{x{) - — ^F^X! -x k ) 

J V k= 1 


N 


p{Xm)da 


m =1 


< 


N 


F * p(xi) — j^'^2F(x 1 -x k ) 


N 

2 p/2 


k =1 


N 


p/2 


p{x m )d 2 


m =1 




Appendix B. Toplitz operators, Wigner and Husimi transforms 

For each z £ C d with 5ft(z) = q and 5s(z) = p, we denote by |z,e) the wave 
function (sometimes referred to as a “coherent state”) defined by the formula 

(45) \z, e){x) := (Tre)-^-^-^ 2 ^'^ . 

We recall the bra-ket notation: |z,e)(z,e| designates the orthogonal projection on 
the line C|z,e) in L 2 (C d ). 

With the normalization above, denoting fj := L 2 (R d ), one has both 
(z,e|z,e) := |||z,e)||; = / | \z, e)(x) / dx = 1, 

and 

(46) 1 d [ \z,e){z,e\dz = I# , 

(27re) d J C d 

where the integral on the left hand side is to be understood in the weak operator 
sense. 

To each positive or finite Borel measure p on C d , we define the Toplitz operator 
at scale e with symbol p by the formula 

OP J(P) : = (0 1 w [ \z,e){z,e\p(dz). 

(27re)“ J C d 

This is a possibly unbounded operator on fj, defined by duality by the formula 
HOP J(p)u)sj ■■= j u(x)OPj{p)v(x)dx 

= (2^p J cd (v\z,e){z,e\u)p(dz) 
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for all u, v G Sj such that z i-A (z, e\u) and z i-A {v\z, e) belong to L 2 (C d , p). 

If p is a positive measure, then 

(47) OPf (p) = OPf (p)* > 0 , trace(OP £ T (p)) = 1 d [ n(dz). 

( 27 re)“ J cd 

If / G L°°(C d ), then 

0 p f(/) = 7^—Td / l*> e ><z> e l/(~) d ~ e C(Sj) , With II OPf (/)|| < ||/|| L =o. 
(2ire) J C d 

The following formulas are elementary but fundamental: if / is a quadratic form 
on R d , then 

f OPf(/(g)) = (/(*) + |e(A./)/*), 

I OPf (/(p)) = (f(~ied x ) + |e(A/)/*), 


where /(#) designates the unbounded operator defined on P 2 (R d ) by the formula 

(f(x)<t>)(x) = f(x)<l){x) . 

Let A be an unbounded operator on P 2 (R d ), and assume that its Schwartz kernel 
k A belongs to <S'(R d x R d ). In other words, A is the linear map from <S(R d ) to 
tS/R^) defined by the formula 


(Au, f)s'(R d ),S(R d ) — (/C4, f ® M)s'(R d xR d ),S(R d xR d ) • 

The Wigner transform of A at scale e is defined as 

(49) W e [A] := (2Tr)~ d J r 2 (k A o J e ), where J e (x, y) = {x + \ey, x - \ey) 

and where Tq is the partial Fourier transform in the second variable. When k A o J £ 
is an integrable function, one has 

W<\A}(x,£) = (27re) _d [ e^' v/e k A (x + \y, x - \y)dy. 

J R d 

In particular, for each q,p G R d , one has 


(50) 


= (27re) _d 



W e [\q + ip, e)(q + ip,e\](x,^) 
^^-d/2 e -i£-y / t e ip-y / e e -(\x+y/2-q\ 2 +\x- V /2-q\ 2 )/2e ^ 

{ne)- d/2 e- lx - q ^ /e {2Tre)- d [ 2/4e dy 

J R d 


= (ne)- d e-^ 2+ ^ 2)/£ , 


since 


-le-Pl 2 /^ 


[ (e/Air) d / 2 e- iv '^- p) e- ev2/A e-\ x - q \ 2/t dy = e" l? 

J R d 

(which is the classical formula for the Fourier transform of a Gaussian density). 
Thus, for each positive or finite Borel measure on R,' / x R d , one has 

(51) w e [OPf (m)] = ^r d G\% * M, 


where G" is the centered Gaussian density on R n with covariance matrix al. With 
p(dqdp) = dqdp, or p(dqdp) = f(q)dqdp, or p(dqdp) = f(p)dqdp , where / is an 
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arbitrary quadratic form on R ci , one finds that 

f W e [Isi\(x,Q = (2ne)- d , 

(52) lw e [f{x)](x,£) = {2ne)- d f(x), 

{w e [f(-ied x )](x,Z) = (2ne)- d f(Z). 

Indeed, when p(dqdp) = g(q,p)dqdp with g a polynomial with degree at most m, 
one has 


(53) 


r 2d 4 . 

U e/2 


— pt&x.t/ 4 


9 = e 


9 = 


Z- _ j 4 n 72! 

0<n<m/2 


rAI 




Thus, for g of degree < 2, one finds that G 2d 2 * Xi ^g = (1 + |eA x ,i)g- Together with 
(46) and (48), this formula justifies (52). 

The Husimi transform of an operator A at scale e is defined in terms of the 
Wigner transform of A by the formula 

W e [A] := * X)€ W e [A]. 

Let R £ U(f>); for each ifi £ L 2 (R d ), one has 


traceflV’XV’l#) = ( 2 ^) d JJ W e [ip\(x, £)W e [R\(x, £)dxd£ 


by the definition (49) of the Wigner transform, and Plancherel’s identity. Special¬ 
izing this formula to if) = | z,e), one finds that 

traced^, e)(z, e\R) = (z, e\R\z, e) 

= (2ire) d W e [R](q,p), where q = SR(^) and p = $$(z). 

More generally, if p is a positive or finite Borel measure on C d , one deduces from 
the previous identity and the Fubini theorem that 

(54) trace(OP^(/z)i?) = [ W e [R](z)p(dz). 

Jc d 

(The previous formula is the particular case where p = (2ne) d S z .) 

In particular, if R = OPj ((2ire) d p) with p £ V- 2 (R d x R d ), for each quadratic 
form f on R d , one has 

trace((/(ar) + f(-ied x )) OPf ((2ne) d p)) 

= [[ (/(?) + f(p))p{dqdp) + feA/ . 

J J R d xR d 

Indeed, applying (54) shows that 

trace((/(a;) + f(-ied x )) OP^(( 2ne) d p )) 

= trace(OP^( f{q) + f{p) - §e(A/)) OPf ((2 tt e) d p)) 

= [[ (G 2 e d 2 *W e [OPj((2Tre) d p)})(q,p){f(q) + f(p) - \e(Af))dqdp 

J J R d xR d 

= [f W € [OpT((2Tre) d p)](q,p)(G 2 e ‘f 2 *(f{q) + f(p)-±e{Af))dqdp 
JJ R d xR d 

= [[ ( G e/2 * G £/2 * (/(«) + fip) ~ 3 £ (A f))P-( d Q d P) , 

JJ R d xR d 
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where the last equality follows from (51). We conclude by observing that 

G\% * * ( f{q ) + f(p) - ie(A/)) = + f(p ) - ±e(A/)) 

= (f(q) + f(p) + H A f )) 

according to (53). 
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